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Abstract 

Several extensions of the standard electroweak model allow new four-fermion interactions V a vpf f 
with strength e a pGp, where (a, (5) are flavor indices. We revisit their effects on flavor oscillations 
of massive (anti) neutrinos in supernovae, in order to achieve, in the region above the protoneutron 
star, an analytical treatment valid for generic values of the neutrino mixing angles (u, (ft, tft) = 
{6 12, $13) $23)- Assuming that e Q/ g « 1, we find that the leading effects on the flavor transitions 
occurring at high {H) and low (L) density along the supernova matter profile can be simply 
embedded through the replacements (ft —>(/> + en and u> — > uj + e£, respectively, where en and €l 
are specific linear combinations of the e^'s. Similar replacements hold for eventual oscillations 
in the Earth matter. Prom a phenomenological point of view, the most relevant consequence is a 
possible uncontrolled bias (<f> — > <fi + en) in the value of the mixing angle (ft inferred by inversion of 
supernova neutrino data. Such a drawback, however, does not preclude the discrimination of the 
neutrino mass spectrum hierarchy (direct or inverse) through supernova neutrino oscillations. 

PACS numbers: 14.60.Pq, 14.60.St, 13.15.+g, 97.60.Bw 
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I. INTRODUCTION 



Atmospheric and solar neutrino data provide convincing evidence for lepton flavor non- 
conservation in the channels — > v T and v e — > v^ T , respectively The Maki-Nakagawa- 
Sakata-Pontecorvo hypothesis of v masses and mixings suffices to explain such phenomena 
in terms of flavor oscillations |IJ, with no need for additional neutrino properties or interac- 
tions. 

However, new neutrino interactions are often predicted (together with neutrino masses 
and mixings) by several extensions of the standard electroweak theory, in the form of effec- 
tive (low-energy) four-fermion operators O a p ~ v a vpf f with strength G^g, inducing either 
flavor-changing (a ^ (3) or flavor non-universal (a = (3) neutrino transitions. Models with 
left-right symmetry || or with supersymmetry with broken .R-parity j|, [5| provide widely 
studied realizations of such operators. Their net effect in ordinary matter (with fermion 
density Nf(x) at position x) is to provide extra v interaction energies || proportional to the 
following dimensionless couplings 

e aP (x) = £ G f a0 N f (x)/G F N e (x) (1) 

f=e,u,d 

= e e ap + e u a3 N u (x)/N e (x) + e d a3 N d (x)/N e (x) (2) 
= e e ap + e% B - 4p + (e^ + 2e d a3 )/Y e (x) , (3) 

where e^g = G^/Gf, and Y e (x) is the electron/nucleon number fraction. 

As a consequence, in the flavor basis v a = (z/ e , z/^, z/ r ), the standard interaction hamilto- 
nian in matter || ^j 1 

H«* = diag(^ 0, 0) , (4) 

can be generalized as 

(H int ) a3 = (H^t)a/3 + Ve a3 , (5) 

where 

V(x) = V2G F N e {x) . (6) 

The interaction hamiltonian for antineutrinos is obtained by replacing G* aB with — (G^g)*, 
namely, V with —V and e a p with e* ag . 2 

Nonstandard neutrino interactions in matter, and their interplay with the oscillation phe- 
nomenon, have been investigated in many different contexts. An incomplete list includes 
analyses related to the solar neutrino problem ||tg[[g§ [13], |TJ, [TJ, [T§, [T7|, |IJ, |TJ [2(| [21 



to the atmospheric neutrino anomaly ||22|| , to the LSND oscillation results [Q, to neutri- 
noless double beta decay [P^j , to supernova neutrinos £31, to the production or de- 
tection of laboratory neutrinos [^, ^7], ^HJ, and to future (very) long baseline projects 



[50 , |3"T] , |3"2| |3B| , P^ , |55| . In a nutshell, these works demonstrate that nonstandard neu- 
trino interactions in matter can (sometimes profoundly) modify the interpretation of current 
and future oscillation searches, and that their possible effects deserve particular attention, 
whenever inferences about standard mass-mixing parameters are made from experimental 
data. 



1 Up to terms proportional to the unit matrix, and discarding negligible one-loop corrections 

2 Notice also that e* aj3 = ep a for hermiticity of H lnt . 
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In this work, we revisit the effects of subleading nonstandard four-fermion interactions 
on the flavor oscillation of massive and mixed (anti)neutrinos in supernovae. We shall limit 
ourselves to the region above the protoneutron star, where neutrinos propagate freely, and 
all relevant fermion densities are monotonically decreasing (roughly as the third power of the 
radius). 3 For this case (previously studied, e.g., in jn| we aim to a completely analytical 
treatment, valid for generic values of the neutrino mixing angles and for both the direct and 
inverse mass spectrum hierarchies allowed by the current v phenomenology. Such treatment 
has been recently proposed in R?| for the case of standard supernova v oscillations. 

In the present work, we make a few simplifying assumptions about the nonstandard 
couplings e a /3- (1) they are taken as real (either positive or negative), so that e a p = ep a ; 
(2) their possible variation with x through Y e [see Eq. (|||)] is neglected; 4 and (3) they are 
assumed to be small, 

M < 1 , (7) 

sa y, \ e ap\ < f ew x 10 -2 , so that terms of 0(e 2 ) can be safely neglected in the calculations. 

Under the above assumptions, we find that the main effect of nonstandard interactions on 
supernova neutrino oscillations is to induce a shift in the relevant mixing angles, as compared 
to the standard case. Such a shift is discussed first in the two-neutrino case (Sec. II) and 
then in the three-neutrino case (Sec. III). Nonstandard effects in the Earth matter are briefly 
discussed in Sec. IV. Phenomenological consequences are illustrated in Sec. V. Conclusions 
are drawn in Sec. VI. 

A final remark is in order. The assumption in Eq. (|7|) is consistent with phenomenological 



upper bounds |L8|, [19|, |3| |40fl , derived from charged lepton processes under SU(2) symmetry 
assumptions. Such bounds are typically of O(10 _1 ) for the (differences of) diagonal couplings 
and of O(10 -2 ) for the off-diagonal ones; only the upper limit on e eAt is much more severe 
[0(1CT 5 )]. However, if allowance is made for substantial SU(2) breaking, and if the analysis 



is performed in a more model-independent way, the previous constraints [18|, |19| |23|, f4C§ can 
be significantly relaxed [41, 42| . Accordingly, scenarios with (some) large e Q/ g's [~ few x 10 _1 ] 



can be built (see, e.g., |41], [43|]), but they will not be considered in this paper. 



II. TWO-NEUTRINO TRANSITIONS 

In this section, we study the effects of nonstandard interactions on 2u oscillations, in 
the region well above the neutrinosphere. We show that the 2u evolution equation can be 
formally cast in an effective (and solvable) standard form (e a! 3 = 0), through a proper redef- 
inition of the neutrino mixing angle. The traditional formulae, derived through resonance- 
condition arguments |L6|, [17], [18], |25| (inapplicable for mixing angles > 7r/4, see |38j and 



references therein), are shown to be recovered in the limit of small mixing. At the end 
of this section, we also comment on a different mechanism that might induce nonstandard 
flavor transitions near the neutrinosphere. 

In the flavor basis, the two-family evolution equation for free-streaming neutrinos can be 



3 Possible flavor transitions induced by nonstandard interactions near the neutrinosphere |3^, [37J will be 
commented at the end of Sec. II. 

4 Detailed supernova simulations show that Y e ~ 0.5 above the protoneutron star. Significant variations 
(reductions) of Y e are confined to regions near the neutrinosphere during the cooling phase 



3 



written as 




[HuaV) + Htoix)] [ ^ ) , (8) 



a 



where a = \x or r, and H = H]^ n {9) + H int (x) is the total hamiltonian, split into kinetic and 
interaction energy terms. The kinetic term can be expressed in terms of the neutrino mixing 
angle 9, of the squared mass difference Am 2 = m\ — m 2 , and of the energy E as 

H km = U(6)(- k Q /2 + l /2 ^jU\6), (9) 

where k = Am 2 /2E is the neutrino oscillation wave number, and 

/ cos 9 sinflA 
w \ — sin 9 cos 9 I y ' 

is the mixing matrix. The interaction term can be written as 

H int (x) = V{x)( l+ c €ee e f a ) (11) 




-1 + t ee + e a 



(12) 



where we have made H [nt traceless in the last equation, V(x) being defined as in Eq. (^|). 
For e a p = (standard case), the interaction term reduces to 

H?£ = diag(+l, -1) ■ V(x)/2 . (13) 

The term H mt in Eq. (|1^) is diagonalized through 

#int = U( V y diag(+l, -1) U( V ) ■ f V(x)/2 , (14) 

where 



/= v /(l + 6 ee - eaa ) 2 + 46 2 a , (15) 
while U(r]) is defined as in Eq. fliPp, but with argument 

7] = - arctan — (16) 

2 1 ~\~ t-ee ^aa 

= e ea + 0(e 2 ), (17) 

where terms of 0(e 2 ) can be neglected under the assumption in Eq. (0). 

Equation (|14|) differs from Eq. ( |13"1) in two respects: (a) the potential V is multiplied by 
a factor /, and (b) there is an additional rotation matrix U(r]) = U(e ea ) + C(e 2 ). The point 
(a) is equivalent to rescale the supernova density by less than a few percent [/ = 1 + 0(e)], 
and does not lead to realistically measurable consequences; in fact, the difference / — 1 
is conceivably smaller than the uncertainties associated to the normalization and shape of 
inferred or simulated supernova density profiles; moreover, the oscillation physics depends 
on V mainly through its logarithmic derivative |38| . The point (b) can instead be important 



and, as we shall see, it leads to a bias in the "effective mixing angle" governing flavor 
transitions along the supernova matter profile. 
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Summarizing the above discussion, the leading nonstandard effect on 2v propagation in 
supernovae is governed by the off-diagonal coupling e ea , and amounts to write ifi nt in the 
form 

tfint(*) * ^ U\e ea ) ( + 1 \ ) U(e ea ) . (18) 
It is then easy to check that, in the new basis 





(19) 

the neutrino evolution equation in matter [Eq. (|)] takes an effective "standard" form, 

[HU0)+H^(x)}(^) , (20) 




where 9 = 9 + e ea . Roughly speaking, nonstandard effects can be "rotated away" from H int 
(which becomes H^f ) and "transferred" to ifkin through the replacement 9^9=9+ e ea . 
Since accurate analytical solutions of Eq. ( |20| ) exist [expressed, e.g., in terms of the survival 
probability P[y e — > u e )], the task is reduced to a proper back-rotation of known results 
(applied to z/ a ) to the physical flavor basis u a . 

To do so, one defines the usual mixing angle in matter 9 m diagonalizing the total hamil- 
tonian in Eq. (PUD, 

„z cos 20 -Wife 

cos20 m = 1 = (21) 

(cos 20- V/k) 2 + sin 2 29 



and the associated basis of effective mass eigenstates in matter v m ^ 

2>a = U(6 m ) Vm,i • (22) 

The formal solution of Eq. (p0|), evolving the state {v a ) from the starting point (x = x ) 
to the detection point (x = Xd), can then be factorized as 



U{d)-T v .T s {e)-U\e° m )U{e ea )[ ^ (23) 





U(9)-T v .T s (9).tf(9° m -e ea )[ P ° ) , (24) 



o 



where, from right to left, U(e ea ) rotates the initial state u a to the v a basis of Eq. ([20|), 
W(9^) rotates v a to the initial basis of mass eigenstates in matter, T$(9) embeds the 
transitions v^a — > P m ,i U P to the supernova surface (where V = and v m ^ = vi), T} 



v 



further propagates the ^'s in vacuum until they arrive at the detector (up to Earth matter 
effects), and a final rotation U{9) brings the v$ states back to the physical flavor basis v a . 



Equation (24), obtained by grouping the first two rotations, represents the formal solution 
of the v evolution equation, in a format analogous to the standard case (i.e., initial rotation, 
propagation in matter, propagation in vacuum, final rotation). 

One can then average out unobservable oscillating terms (mainly associated to TV), so 
as to propagate classical (real) probabilities P a p rather than complex amplitudes. The 



averaging can be accomplished by replacing in Eq. fl24|) each entry in the matrices 
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U(9), Ts(0), and 1^(9^ — e ea ), with its squared modulus. In our case, the initial (high 
density) condition V(x Q )/k 3> 1 leads to sin 2 (6^ — e ea ) ~ 1 + 0(e 2 a ), effectively cancelling 
any "memory" of the initial point x D above the protoneutron star. Therefore, the electron 
neutrino survival probability finally reads 5 

_ r x / cos 2 fl sin 2 # \( 1 - P c {9) P c {9) WoiWlN 
ee ~ V 1 U J ^ sin 2 9 cos 2 9 ){ P c (§) 1 - P c (9) ) { 1 ) { ) ' ^ 

namely, 

P ee = cos 2 9P C {9) + sin 2 0[1 - P c {9)) , (26) 

where we have adopted the usual notation for the square moduli in T s , in terms of the 
so-called crossing probability P c (see, e.g., J44|). 

Equation (|26|) is formally analogous to the one in the standard case (e ea = 0), modulo the 
replacement 9 — > 9 = 9 + e ea in P c . We recall that, in the standard case, P c {9) is accurately 
described by a simple analytical expression, 

pl-Krk cos 2 9 i 

= e „r k _ I » (27) 



where 

r = — 



dlnF 

dx 



(28) 



is the density scale factor, to be evaluated at the point x p where the potential equals the 
wave number |38|] , 

V(x p ) = k . (29) 
The above prescription for P c std (#), inspired by the condition of maximum violation of adi- 



abaticity [p8| , |47[] (more general than the traditional condition of resonance ^3|) and by 



the double-exponential representation for the crossing probability | i8fl , holds accurately for 
generic supernova density profiles, in the whole 2u mass-mixing parameter space p8| . 

In the presence of nonstandard interactions, the previously discussed replacement 9 — > 
9 + e ea leads then to the following analytical form for P c : 



2irrkcos 2 (9+e ea ) _ i 

PM = ^ _ i — . (so) 



where r takes the same value as in Eq. (p8|), the position of x p [Eq. (p9|)1 being unaffected 
by the shift 9^9 + e ea . 6 

Equations (|26|) and (|30|) represent our main result for the two-flavor neutrino case. They 
provide a compact description of nonstandard interaction effects on P ee , valid for generic 
values of 9 and Am 2 (provided that e ea 1). The calculation of P ee for antineutrinos 



5 We remind the reader that P ee is the most relevant quantity in supernova v oscillations, due to the 
practical indistinguishability of and v T in the initial and final states [Q. See, however, Ref. |^5| for a 
recent discussion of possible differences between and v T fluxes at the origin. 

6 The definition of r(x p ) through Eqs. ( ^8| ) and (|2^ ) is ^-independent. 
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is strictly analogous: given the analytical expression reported in Sec. II. C of |38|] , one has 
simply to replace 9 with 9 + e ea within P c (u). 7 

It can easily be shown that Eq. ( |30"D generalizes previous results about nonstandard 
interaction effects, valid in a more limited range of applicability. The traditional analytical 
approach to P c localizes the nonadiabatic transitions at the so-called "resonance" point 
x r , where mixing in matter is maximal (sm29 m (x r ) = 1) |TB|, [T7|, |TB, 2E\. The resonance 



approach is certainly valid at small 9, where it should coincide with ours. Indeed, in the 
resonance approach, the standard (e ea = 0) Landau-Zener form for P c at small mixing 
(lnP c oc 9 2 ) g§ is modified as lnP c oc 9 2 {l + 2e ea cot2#) 2 ~ (9 + e ea ) 2 |]| |T7|, [18], in 
agreement with our prescription 9 — > 9 + e ea . For 9 approaching 7r/4, however, the resonance 
condition generally fails to localize the right point x for the evaluation of r(x), and becomes 



eventually not applicable for 9 > n/A (both in the standard |38] and in the nonstandard 
case). Conversely, Eq. ( |30D has no such applicability restrictions. Finally, in the subcase 
9 = 0, we also recover from Eq. (BO) the double-exponential form for P c found analytically in 



|TT|| in the scenario of massless (solar) neutrinos with nonstandard interactions. Therefore, 



Eq. ( p0|) correctly generalizes previous limiting cases. 

A final comment is in order. Additional nonstandard flavor transitions (different from 
those described in this section) might independently occur in inner layers 0, |37j| , as a 
result of a decrease of Y e (x) to 0(1CT 2 ), localized just above the neutrinosphere during the 
cooling phase |37|, |39f . The e Q/ g(x)'s can then locally increase to 0(1) [see Eq. ([[])], making 



the diagonal and off-diagonal terms in the interaction hamiltonian [Eq. (|TT|)1 comparable 
in magnitude, and dominant (due to the very high matter density) with respect to the 
kinetic terms [Eq. (61)]. As a consequence, near the neutrinosphere (and thus prior to the 
flavor oscillations discussed in this paper), nonstandard interactions might induce additional 



flavor transitions of (effectively massless) v and v during the cooling phase |36|, |3"? |. Our 
results, applicable to the region of decreasing fermion density above the protoneutron star 
(where Y e ~ 0.5) should be thus considered as complementary to the possible independent 
effects [F7| possibly induced by a local (preceding) drop of Y e . A combination of the two 
nonstandard interaction effects (still lacking in the literature, to our knowledge) is, however, 
beyond the scope of this paper. 8 



III. THREE-NEUTRINO TRANSITIONS 

In the 3z/ case, the interaction term in the Hamiltonian is given by Eq. (^), while the 
kinetic term H^ Q is 

f-k L /2 \ 
= U(u, <f>, V) +k L /2 \u\u,4>,^), (31) 
V kfj J 



7 Notice that, outside P c , the angle 8 is unshifted, as in the coefficients sin 2 8 and cos 2 8 of Eq. (|26|). 

8 We note in passing that, being confined to a region just above the neutrinosphere (where collective and 
collision effects can still be non negligible), the additional transitions at low Y e might require a more 
advanced description than the free-particle formalism used in []37| (e.g., in terms of kinetic and density 
matrix equations |l9| ) . 
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where ki = Sm 2 /2E and kn = m 2 /2E (6m 2 and m 2 being the "solar" and "atmospheric" 
neutrino squared mass differences, see |38|), while the 3x3 mixing matrix U(u, <f), if>) is 
factorized (discarding a possible CP violating phase) into three real rotations H, |5(J 

U(u, 0, V) = U^)UM)U n {u) , (32) 

where (u,<p,i/j) = (#12, #13, $23)- The subscripts H and L remind that two transitions can 
occur along the supernova profile when V ~ O(kn) and V ~ 0(ki,), namely, at high (H) 
and low (L) density, respectively. 

The phenomenological assumption of hierarchical squared mass differences, 

5m 2 < m 2 < > k L <^k H , (33) 

together with Eq. (0), makes the 3z/ dynamics factorizable into two (almost decoupled) 2v 
subsystems for the H and L transitions, in the same way as for the standard case H . 

To isolate the dynamics of the H transition, one rotates the starting neutrino (flavor) 
basis by U^ip), and extracts the submatrix with indices (1, 3) [ER 53]. The effective 2v 



sub-dynamics is then governed by (k, 9) = (kn, 0), and the analogous of e ea is given by the 
off-diagonal (nonstandard) term, which reads 

£h = ten sin ip + e er cos ip . (34) 



Using the results of the previous section and those of Ref . , the crossing probability for 
the H transition is thus described in terms of an effective (1, 3) mixing angle <ft = <ft + e#, so 
that 

e 2nr H k H cos 2 ((p+e H ) _ j 
PH = e 2nr H k H _ I ' ( 35 ) 

th being the density scale factor [Eq. (p5[) ] evaluated at the point xh where V(xh) = ku- 

Analogously, to isolate the dynamics of the L transition, one rotates the starting neutrino 
(flavor) basis by £^13 (0) C/23 (V*) 5 an d then extracts the submatrix with indices (1, 2). The 
effective 2u sub-dynamics is then governed by (k, 8) = [ki,,uj) and by the off-diagonal 
(nonstandard) term, which now reads 9 

ei = cos (fi(e efJi cos ip — e er sm ip) 

+ sm(f)[(e TT — e w ) sin ^ cos ^ — e (UT cos2'0] . (36) 

The crossing probability for the L transition is thus described in terms of an effective (1, 2) 
mixing angle £0 = uo + 6l, so that 

e 2irr L k L cos 2 (w+e L ) _ -y 
PL = e 2*r L k L _ I ' ( 3? ) 

r L being evaluated at the point xl where V(xl) = k L . The same shift (uj — >• u + e^) applies 
to the expression of Pi for antineutrinos (as given in f38|). 

The above discussion implicitly assumes direct hierarchy of neutrino masses (m 12 < 
777.3). The case of inverse hierarchy (m 3 < m 12 ) is, however, strictly analogous; indeed, the 



An analogous result for has been recently reported in the context of solar neutrinos plj ] . The expression 
of el for supernova neutrinos, in the subcase e Mr = e w = e TT = 0, is reported in [ p5| . 
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symmetry arguments used in |38| to derive P ee analytically in the case of inverse hierarchy (in 
terms of P ee for direct hierarchy) are unaffected by the presence of nonstandard interactions. 
Therefore, the prescription remains the same, irrespective of the spectrum type: take the 



standard expressions of P ee from |]3~8fl , and replace uo — > uj + €l and — > + e# in the crossing 
probabilities for both neutrinos and antineutrinos. Embedding nonstandard interaction 
effects through appropriate shifts of u and within Pi and Ph, respectively, represents our 
main results for the case of three-flavor mixing of supernova neutrinos (up to Earth matter 
effects) . 



IV. EARTH MATTER EFFECTS 

In the standard case (e a p = 0), it is well known that possible Earth matter effects 
before supernova v detection can be embedded through a specific transition probability 
Pe = P{yi — > v e ) (see |3l| |5"T] , |52 | and references therein). In the hierarchical approximation 
[Eq. (|5BD], P E is basically a function of the //-transition parameters (u,5m 2 ) only. 10 

The dependence of P E on u is both explicit (through an tu-rotation in vacuum), and 
implicit [through the mixing angle in the Earth matter, u m = iv m (u)], so that 

Pe — Pe(&, u m (u), 5m 2 ), standard case . (38) 



In particular, Sec. IV of [38] reports known representations of the function 
P e (uj, u) m (u), 5m 2 ) in the case of Earth mantle (+ core) crossing, in the same notation 
as in the present paper. 

In the nonstandard case, since it is e a p ^ only in matter, the previously found shift 
uj — > uj + €l applies only to u m , namely, 

Pe = Pe(u, co m (tu + €l), 5m 2 ), nonstandard case , (39) 

and similarly for antineutrinos. 11 We have also checked the above prescription by explicitly 
repeating the derivation of P E in the cases of one-layer or two-layer model for the Earth 
crossing (omitted). 



V. DISCUSSION 

Let us discuss in more detail the effects of nonzero and en on supernova neutrino 
oscillations, by considering the phenomenologically interesting case with m 2 = 3 x 10~ 3 eV 2 , 
5m 2 m 2 , and sin 2 < few% |1[|. We assume the supernova power-law density profile 
reported in Fig. 1 of |3SJ, and a representative neutrino energy E = 15 MeV. The electron 
neutrino survival probability, in the case of direct hierarchy, is then given by 

P ee {v) ~ [(1 - P E ) P L (k L ,uj + e L ) +P E {1- P L (k L ,u + e L ))\ e -**rx*B ? ( 40 ) 



There is a subleading dependence on through the renormalization V —> V cos 2 <f>, which can be safely 
neglected for supernova neutrinos, given the smallness of <fi Q and unavoidable normalization uncertainties 
of V. 

In principle, the value of cl to be used in Eq. ( p^ ) might be slightly different from the one in Eq. (|37|), 
due to different (average) values of Y e in the Earth and in the supernova [see Eq. (||)] . 
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where we have neglected a small (additive) 0(4> 2 ) term, the main effect of being embedded 
in the exponential suppression factor. In the above equation, Pl is given by Eq. (|37|), while 
Pe is calculated as in Eqs. ( |38]) or (|39|). In the case of no Earth crossing, Pe is simply given 
by 

Pe = sin 2 co> (no Earth effects) . (41) 

Figure 1 shows curves of iso-P ee (z/), calculated according to Eq. (0), with and without 
Earth matter effects (dotted and solid lines, respectively), 12 in the plane of the L-transition 
variables (5m 2 , tan 2 u). The three left panels refer to the case tan 2 (0 + e#) = (correspond- 
ing to either = en = or to tu = ~<P), which implies Ph = 1. The three right panels 
correspond instead to a nonzero value for + e#, namely, tan 2 (0 + e#) = 2 x 10 -5 , implying 
Ph = 0.46. The comparison of left and right panels clearly shows the suppression of P ee {v) 
due to + tft 7^ 0. 

In Fig. 1, from top to bottom, the three (left and right) panels correspond to = 0, 
6l = +3 x 10~ 2 , and = —3 x 10~ 2 . The effect of nonzero is dramatic at small 
mixing angles, and leads to a distinctive and well-known pattern for the isolines of P ee 
T2| , |T4| , |TT| , |T5| , |25|1 , which appear to be elongated towards very small values of uj for positive 



€l, and split into two branches for negative e^. In our formalism, this behavior is simply 
captured in terms of the cos 2 (u; + e^) argument in Pl. for > the argument increases 
monotonically for decreasing (small) uj, otherwise it first reaches a maximum at uj = — 
and then decreases. The Earth effect is also similarly "stretched" in Fig. 1. 

Let us consider in more detail the case of relatively large uj [tan 2 uj ~ O(10 ±:L )], which ap- 
pears to be generally favored by recent solar neutrino data |53| . In this case, independently 
of the value of 5m 2 , effects induced by ^ are hardly recognizable in Fig. 1, since the 
difference between uj and uj + is fractionally small, while the leading effect of nonstandard 
interactions is embedded in the exponential suppression factor (Ph) of Eq. ([4(]). In partic- 
ular, for relatively large values of 5m 2 / E [providing the best fit to solar neutrino data |53| 



through the so called large mixing angle (LMA) solution], it is Pl ~ 0, and the v e survival 
probability simply reads 

P ee ~ P E e -2^„fcffsin 2 «,+e H ) (LMA solution) . (42) 

Figure 2 shows the function P E (5m 2 ) for three representative values of tan 2 ci; in the 
LMA region. Other parameters (e^, neutrino energy and pathlength, mantle density and 
electron number fraction) are as in Fig. 1. The case of no Earth crossing [Eq. (|4"T|)1 is 
also shown as a horizontal line in each panel. In general, Earth matter effects appear to 
be relevant in all cases, producing a distinctive oscillatory pattern in P E . Different values 
of €l induce slightly different oscillation amplitudes in P E , and could thus (in principle) 
be distinguished, especially for large values of 5m 2 , where the relative differences among 
the curves in each panel increase. However, such discrimination would require very high 
statistics observations. We conclude that, for mass-mixing parameters in the LMA solution, 
and for typical supernova density profiles and energy spectra, the main nonstandard effects 
are associated with en (through Ph), while smaller effects are associated with (through 
Pe). 



12 In Fig. 1, Earth matter effects refer to a representative case of mantle crossing, with pathlength L = 8500 
km, mass density p = 4.5 g/cm 3 , and electron number fraction Y e = 1/2. 
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Concerning the main nonstandard effects, the survival probability in Eq. ([40|) shows that, 
in general, supernova neutrinos can be sensitive to values of tan 2 (0 + e#) as low as ~ 10~ 5 , 
corresponding to a (sub)percent sensitivity in + e#. Indeed, in the standard case (e# = 0), 
it has been shown that inversion of future galactic supernova data might provide upper 
and lower bounds on in the percent range ||54j| . Our results show that, in the presence 
of additional nonstandard interactions (e# 7^ 0), any constraint about must actually be 
interpreted as a constraint on + e#, namely, there is a strict degeneracy between the 
vacuum mixing angle = #13 and the nonstandard coupling e# given in Eq. (0). 

The (0, €h) degeneracy in supernova neutrino oscillations can be rather dangerous for 
the interpretation of experimental data. Supernova v oscillations may offer one of the few 
opportunities to probe experimentally very small values of 0, and any uncontrolled bias of 
the kind — > + e# may lead to dramatic differences in the theoretical inferences from 
the experimental data. In extreme cases, nonstandard interactions might either completely 
mimic (0 = 0, e# 7^ 0) or completely cancel (0 = —e H ) "standard" 0-related effects in 
the H transition. 13 An analogous "confusion scenario," with degeneracy between effects of 



nonstandard interactions and of (f> 7^ 0, has recently been discussed in |}5] in the context 
of oscillation searches at neutrino factories. The main message emerging from such results 
is that constraining will generally be a very challenging task, if allowance is made for 
nonstandard (flavor changing) neutrino interactions, with strength as weak as a few percent 
of the standard (flavor diagonal) electroweak one. 14 

We conclude this section with a positive remark, by showing that not everything is 
necessarily biased by nonstandard interactions. In particular, the spectrum hierarchy dis- 
crimination (direct vs inverse) may still be viable. Let us consider the phenomenologically 
interesting case of the so-called large mixing angle solution (LMA) to the solar neutrino 
problem |k|, requiring 5m 2 ~ 0(1CT 5 ) eV 2 and tan 2 u; < 1, so that Pl — (adiabatic 
L transition). For small 0, the supernova electron (anti) neutrino survival probabilities are 
approximately given (up to Earth matter effects) by 

P^(u) ~ sin 2 ue -2nr H k H si^(ne H ) ^ ( 43 ) 

P e d e ir (z7) ~ cos 2 cu , (44) 
P*>) ^ shi 2 ^, (45) 
P^ v (z/) ~ cos 2 ue~ 2wrHkHSin2 ^ +eH) , (46) 

where the superscripts denote the cases of direct or inverse hierarchy. From such equations 
one can derive the following e#-mdependent inequalities, 

Pf:^) < , (47) 

Pfe{y) > F^iy) , (48) 

and 



13 



We remind the reader that, even assuming restrictive bounds on the e a p |T^, [19| 
currently exclude values of e# as large as a few percent. 



40 



one cannot 



At a more sophisticated level, for complex e a p there can also be confusion between standard and nonstan- 
dard CP- violation phase effects in (very) long baseline experiments; see, e.g., pl 



:V2 
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which become all stronger for decreasing neutrino energy. Such relations can then be used 
to infer information about the mass spectrum hierarchy from general features of supernova 
(anti)neutrino event spectra, 15 analogously to the standard case |52], [55], [5(J. The hierarchy 
discrimination becomes impossible, however, for sin 2 (0 + ejj) <C 10 -5 (corresponding to 
exp[— 27rr H kH sin 2 (0 + e#)] — 1), in which case the above three relations become equalities. 

We finally remind that, under certain conditions (see the end of Sec. II), additional 
nonstandard transitions might take place just above the neutrinosphere ff7[] . In this 
case (not considered here), the global P ee function should be obtained by convolving such 
transitions with those (H and L) considered in this and in the previous section, leading to 
a more complicated phenomenology. 

VI. SUMMARY 

We have revisited the effects of nonstandard four-fermion interactions (with strength 
^^Gf) on supernova neutrino oscillations occurring above the protoneutron star. Under 
reasonable approximations, we have found that the oscillation probability can be written 
in an analytical form, valid in the whole 2v parameter space, and applicable also to the Zv 
parameter space in hierarchical (direct or inverse) cases. Our approach generalizes previous 
results, which are recovered as specific limits. We find that, as far as the transitions at 
high (H) and low (L) density are concerned, the main effects of the new interactions can 
be embedded through (positive or negative) shifts of the relevant mixing angles <ft — $13 
and uj = 9i2, namely, — > <fi + en and u — > u + €l, respectively [see Eqs. (0) and fl36f) 1. 
Barring the case of small u (disfavored by solar neutrino data), and apart from small el- 
induced Earth matter effects, the main phenomenological implication of such results is a 
strict degeneracy between standard (0) and nonstandard (e#) effects on the H supernova v 
transition. However, such a degeneracy does not necessarily spoil the discrimination between 
direct and inverse neutrino mass spectrum hierarchy. Our work is complementary to studies 
|37j of nonstandard transition phenomena possibly occurring just above the neutrinosphere. 
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3v surviva I probability, E — 1 5 MeV, m 2 - 3 x 1 0" 3 eV 2 
tan 2 (^ + £ H ) = tan 2 (<^ + £ H ) = 2 x 1 0" 5 




10 10 10 10 1 10 10 10 10 10 10 10 1 10 
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FIG. 1: Isolines of P ee for supernova neutrinos with direct mass hierarchy, in the plane of the 
L-transition variables {5m 2 , tan 2 u), with (dotted lines) and without (solid lines) Earth effects 
(calculated for L = 8500 km, p = 4.5 g/cm 3 , and Y e = 1/2). Nonstandard interactions are 
parameterized through the dimensionless couplings 6l,h- The values of m 2 and E are fixed at 
3 x 10~ 3 eV 2 and 15 MeV, respectively. Left (right) panels refer to tan 2 (0 + e#)=0 (= 2 x 10~ 5 ). 
The parameter ei takes the representative values (upper panels), +3 x 10~ 2 (middle panels), and 
—3 x 10~ 2 (lower panels). See the text for details. 
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FIG. 2: Earth matter effects: Curves of the probability Pg = P{^2 — * v e) as a function of 5m 2 
for mantle crossing of supernova neutrinos (with E, L, p and Y e as in Fig. 1). The chosen range 
of 5m 2 and the three representative values of tan 2 uj refer to the so-called LMA solution to the 
solar neutrino problem [^]. Nonstandard interactions affect the amplitude and the phase of the 
oscillatory pattern through e^. The horizontal solid line represents the case of no Earth effect. 
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